Efficiency of initiating cell adhesion in hydrodynamic flow 
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We theoretically investigate the efficiency of initial binding between a receptor-coated sphere 
and a ligand-coated wall in linear shear flow. The mean first passage time for binding decreases 
monotonically with increasing shear rate. Above a saturation threshold of the order of a few 100 
receptor patches, the binding efficiency is enhanced only weakly by increasing their number and 
size, but strongly by increasing their height. This explains why white blood cells in the blood flow 
adhere through receptor patches localized to the tips of microvilli, and why malaria-infected red 
blood cells form elevated receptor patches {knobs). 



Cohesion in biological systems and biotechnological ap- 
plications is usually provided by specific bonds between 
receptors and ligands. The formation of these bonds re- 
quires a physical transport process which brings recep- 
tors and ligands to sufficient proximity for binding. On 
the cellular level, one of the most prominent examples is 
the binding of blood-born cells to the vessel walls under 
the conditions of hydrodynamic fiow. For white blood 
cells, initial binding to the vessel walls is the first step in 
their hunt for pathogens, which is then followed by rolling 
adhesion, firm arrest and extravasation jlj. Similar pro- 
cesses are used by stem and cancer cells which travel the 
body with the blood stream. Initiating binding to ves- 
sel walls is also essential for malaria-infected red blood 
cells in order to avoid clearance by the spleen and pos- 
sibly also to foster rupture and release of new parasites 
into the blood stream _2J . Similar questions about initial 
binding under flow conditions arise for bacteria, e.g. when 
binding to the intestinal wall ^ , and in biotechnological 
applications, e.g. for adhesion-based cell sorting 4J. In 
order to control shear fiow and cell density in a quanti- 
tative way, a standard setup are flow chambers An 
essential but largely unexplored aspect of these processes 
is receptor geometry, that is size, height and separation 
distance of the receptor patches. One way to address this 
issue experimentally is the use of receptor-coated beads 



Earlier theoretical efforts in this context have been fo- 
cused mainly on issues related to white blood cells, in- 
cluding modelling of the initiation of adhesion at high 
cell densities (e.g. due to hydrodynamic interactions) 7j 
and the process of rolling adhesion |8|. In these stud- 
ies, the parameters characterizing receptor-ligand bind- 
ing are usually fixed at physiologically motivated values. 
In this Letter, we take a more general view and ask how 
a variable receptor geometry affects cell capture in hy- 
drodynamic ffow. In order to study this issue in a sys- 
tematic way, we spatially resolve receptors and ligands. 
The cell is modeled as a rigid spherical Brownian par- 
ticle in linear shear flow carrying receptors for ligands 
covering a planar boundary wall. In the absence of inter- 
actions with other particles or external forces, there is no 



FIG. 1: A cell of radius R moves in linear shear flow above a 
wall at a height z, subject to hydrodynamic, gravitational and 
thermal forces. The efficiency for initiating binding is assessed 
by calculating the mean first passage time for approach of 
receptor and ligand to a capture distance ro, that is for overlap 
of the gray regions. Three models of increasing complexity 
are studied here: (a) homogeneous coverage of cell and wall, 
(b) Nr equidistantly placed receptor patches on the cell, each 
with radius r^, and (c) ligand patches of radius separated 
by distance d. 



reason for such a particle to drift towards the wall and 
initial binding has to rely completely on thermal diffu- 
sion. In order to arrive at a generic model, we consider 
the simplest of the possible downward driving forces in 
experiments with flow chambers, namely gravity. The 
particle is set free at a certain height above the wall and 
we calculate the mean first passage time (MFPT) for the 
first receptor-ligand encounter as a measure for the effl- 
cieny of initial cell binding. We consider three models of 
increasing complexity in regard to the spatial distribu- 
tion of receptors and ligands. We first show that if the 
receptors on the cell and the ligands on the substrate are 
distributed homogeneously, then the corresponding MF- 
PTs can be calculated exactly. In the case that receptors 
are spatially resolved, we use extensive computer sim- 
ulations to calculate the MFPTs as a function of their 
number and spatial dimensions, both in two and three 
dimensions. As a third case, we in addition consider spa- 
tially resolved ligand distribution. 

Fig.nintroduces the parameters of our model. We con- 
sider a sphere of radius R which moves with the hydro- 
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dynamic flow in positive x-direction at a height z above 
the wall with normal e^. The simplest possible flow pat- 
tern is linear shear flow with shear rate 7. With the 
no-slip boundary condition at the wall, the unperturbed 
velocity profile reads u°° = '^ze^. For a typical cell ra- 
dius R = 5 fim and a typical shear rate 7 — 100 Hz, 
the Reynolds number in aqueous solution is well below 1 
and the hydrodynamic flow is essentially described by the 
Stokes equation for incompressible fluids. This is even 
more true for smaller particles like micron-sized beads. 
Scaling estimates also show that for typical parameter 
values for cell elasticity, deformations due to shear flow 
and lubrication forces are small and therefore the spher- 
ical approximation is justified. In addition to the hydro- 
dynamic forces, in our model there are also gravitational 
and thermal forces acting on the particle. Such a combi- 
nation of forces is the subject of Stokesian Dynamics 
and in our case leads to the following Langevin equation 
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FIG. 2; Effect of shear rate on mean first passage time T 
for initial binding for spatially resolved receptors in (a) two 
(2D) and (b) three dimensions (3D). The scales for length 
and time are R and Q-rtiR? /ksTa, respectively. The Peclet 
number Pe = dnvqR^'y /ksTa and the Peclet number in z- 
direction Pe^ = ing ApR"^ / {SksTa) represent the strengths 
of the hydrodynamic and gravitational forces, respectively, 
relative to the thermal forces, zo = 2, Pcz = 50, ro = Vp = 
10"^ 



(1) 



Ta is ambient temperature. The gravitational force reads 
F*^ = — gAp(47ri?'^/3)e3 with some density difference Ap. 
The subscript t denotes random variables. As usual, the 
thermal force gt is assumed to be Gaussian: 







{gtgf) ^ 2kBTaMS{t ~ t') . (2) 



The superscript S for the multiplicative noise term in 
Eq. (QJ indicates that it has to be interpreted in the usual 
Stratonovich sense. The matrix B in Eq. is related 
to the mobility matrix M through M = BB"^. Vectors in 
Eq. 1^ are six-dimensional, representing the spatial and 
orientational degrees of freedom. The mobility matrix M 
and the shear force F'^ for a spherical particle above a 
wall cannot be obtained in analytically closed form. How- 
ever, they can be calculated numerically to high accuracy 
and we will use this for our simulations 

Considering the physical dimensions of our problem 
shows that the motion of the cell is essentially governed 
by two dimensionless numbers. For length, the nat- 
ural scale is cell radius R. Mobility and shear force 
scale as M = l/{6TTr]R) and = GnriR^j, respec- 
tively. For time, there are two relevant time scales, 
the deterministic time scale I/7 and the diffusive time 
scale R^/D = 6nT]R^ /kBTa, where we have used the 
Einstein relation D = MhsTa for the diffusion con- 
stant D. Therefore the relative importance of hydro- 
dynamic to thermal motion is described by the Peclet 
number Pe = %'nr\R?') jkBTa- In the limits Pe and 
Pe 00, diffusive and deterministic motion dominate, 
respectively. The gravitational force introduces another 
dimensionless number, which we call the Peclet number 
in z-direction, Pe^ = F'-^R/ksTa — AirgApR'^ / (3kBTa). 
In the following, we will non-dimensionalize length and 
time by R and 6Trr]R^ /ksTa, respectively. 



We start by considering homogeneous coverage of cell 
and wall with receptors and ligands, respectively, com- 
pare Fig. Then rotational degrees of freedom are 
irrelevant. Because the wall breaks the symmetry only 
in the z-direction, motion in the x-y-plane is decoupled 
from our problem. Thus in this case we essentially deal 
with a MFPT in one dimension, which is independent 
of shear rate 7 and which can be approached with stan- 
dard methods for the appropriate Fokker-Planck equa- 
tion. Binding is identified with approach of receptor and 
ligand to a capture distance tq. Applied to the case of 
homogeneous coverage, the cell has to fall to the capture 
height 1 + rQ. If dropped from the initial height zq, the 
respective MFPT can be shown to be 
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Thus the MFPT scales inversely with the gravitational 
force driving the cell onto the wall. With the lubrication 
approximation Mzz{z) w 1 — 1/z we find 
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Thus the MFPT diverges logarithmically with vanishing 
capture distance ro (that is when the cell has to get in- 
finitely close to the wall) and linearly with initial height 
zq (that is when the cell starts infinitely far away from 
the wall). Although only the constant force chosen here 
results in an analytical result like Eq. Q), for other types 
of force laws it is straightforward to numerically calculate 
corresponding falling times Th [loj . 

We next consider the case of a spatially resolved recep- 
tor distribution, compare Fig.^^- Now the cell is equidis- 
tantly covered with Nr receptor patches, each with radius 
Tp and height ro. We first note that in this case, ini- 
tial orientation becomes important. Moreover now shear 
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rate 7 enters the analysis: the shear flow increases cell 
rotation, and for heterogeneous receptor coverage, this 
strongly influences when the first receptor can bind the 
first ligand. Because experimentally it is hardly possible 
to prepare the initial orientation of the cell, in the follow- 
ing we average over all possible initial orientations. One 
can show that the angle-averaged MFPT is the MFPT 
to fall from initial height zq to some intermediate height 
Zm according to Eq. JSJ (that is independent of orien- 
tation) plus the angle-averaged MFPT to bind from the 
initial height Zm- In this sense, the initial height is not 
relevant for our problem and in the following we always 
use Zq = 2, that is the cell has to fall for the distance of 
one radius before binding can occur. 

In Fig. |2t and b we show the MFPT as obtained by 
extensive computer simulations as a function of Peclet 
number Pe ~ 7 and receptor patch number Nr for two 
(2D) and three dimensions (3D), respectively. Here 2D 
mean that translational motion is restricted to the x-z- 
plane and rotations are restricted about the y-axis, which 
allows for much faster simulations. Each data point in 
Fig.[21is the average of at least 10^ simulated trajectories 
of the Langevin equation Eq. Our simulations are 
very time-consuming because with the receptor patches 
we resolve objects of typical size 10"^, that is nm-sized 
patches on micron-sized cells. From Fig. |21we first note 
that T decreases monotoneously with increasing Pe and 
that the shear rate does not change the relative sequence 
of the curves for different Nr. Thus the larger shear rate 
and the more receptor patches present, the more efficient 
cell capture. The crossover between the diffusion- and 
convection-dominated regimes does not occur at Pe « 1 , 
but at much larger values Pe ~ 10^. Next we note that in 
the 2D case (Fig.|5^), for large Pe or large Nr, all curves 
level off to the exact result for homogeneous coverage 
from Eq. because in these two limits, the binding 
process effectively becomes rotationally invariant. In the 
3D-case (Fig. Eb), the homogeneous reference value is 
only achieved for large receptor numbers. The reason is 
that for small numbers of receptors, the cell might have 
to rotate around the x-axis before a receptor moving on 
a circle parallel to the x-z-plane is able to bind a ligand 
on the wall. Therefore in 3D thermal diffusion remains 
essential even in the case of large Peclet number. 

In order to achieve a better understanding of the simu- 
lation results shown in Fig. 13 it is instructive to decom- 
pose the process into periods of falling and rotation, re- 
spectively. A detailed analysis shows that in the 2D-case 
this decomposition allows to derive scaling laws for dif- 
ferent limits in regard to Pe and Pe^ . An important case 
is the one of large Pe^, when the cell is strongly driven 
onto the wall. Then the binding process can be decom- 
posed into a initial faffing period described by Eq. (PJ, 
followed by a purely rotational search for ligand, which 
is independent of Pe^ and can be calculated analytically: 
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FIG. 3: Mean first passage time T in 3D in the diffusive 
limit (Pe ~ 0) as a function of (a) the number of receptor 
patches A'^^ for different values of patch height ro and (b) 
ligand density pi as varied by decreasing ligand patch distance 
d for different values of the number of receptors iV^. The 
dashed lines show the scaling behaviours T ~ 1/Nr and T ~ 
1/pi at low receptor and ligand coverage, respectively. 20 = 2, 
Pe, = 50, ro = rp = rd = 10"^ 
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Here A9 « 9s is the angle between the absorbing bound- 
aries, 9s = 2iT/Nr the angle between receptor patches, 
Ag K, Pe/2 the rotational drift and Dg « 3/4 the rota- 
tional diffusion constant. From Eq. (0 we get T ^ 1/-/V^ 
and T ~ l/{NrPe) for small and large Pe, respectively, 
in excellent agreement with the scaling found in our sim- 
ulations. In general, Eq. (O is a good qualitative descrip- 
tion of the 2D data shown in Fig. |21 

In order to consider the case of spatially resolved lig- 
and, compare Fig.^, we cover the boundary wall with a 
square lattice of circular ligand patches, with lattice con- 
stant d and patch radius r^. Because again the Peclet 
number does not change the relative sequence of the dif- 
ferent MFPT curves, the efficiency of initial cell binding 
as a function of receptor and ligand geometry can be in- 
vestigated in the diffusive limit Pe « 0. Fig. |21 shows for 
this case that the MFPT saturates both when increasing 
receptor coverage by increasing Nr or ligand coverage by 
decreasing d. A similar saturation behaviour is also found 
when increasing receptor and ligand patch sizes and 
rd- Typically the saturation threshold for the parame- 
ters used is located at a mean patch-to-patch distance d 
of about 0.17, both in regard with receptors and ligands. 
For a coverage below the threshold, Fig.Olshows that the 
MFPT scales like ~ l/TV^ - and - - in re- 
gard to receptor and ligand coverage, respectively, where 
d represents the distance between receptor and ligand 
patches, respectively. This can be understood by noting 
that the ID MFPT for capture by diffusion scales ~ d^ 
where d is the distance between the two absorbing bound- 
aries. The saturation effect observed with respect to Nr, 
pi, Tp and rd results from the space-filling nature of dif- 
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FIG. 4: Mean first passage time T in 3D in the diffusive 
limit (Pe « 0) as a function of tlie receptor patcli radius 
Tp for two different values of ro and three different values of 
Nr. For better comparison with regard to ro the mean first 
passage time Th for homogeneous coverage from Eq. (|3J has 
been substracted. The lines are fits according to Eq. JSJ. 
Other parameters as in Fig.|21 



fusion "which has been imphcated before for the efficiency 
of ligand capture by a cell . 

We finally turn to the effect of the capture distance tq. 
In Fig. 0] "we sho"w the MFPT in the diffusive limit as a 
function of the receptor patch radius rp for ro = 10"'^ 
and rg = 10~^ as "well as for three different values of Nr- 
We first note that the MFPT is much more influenced by 
a change in ro or Nr then by a change in Vp. One can 
sho"w on geometrical grounds that the t"wo parameters rp 
and Tp conspire to define an effective receptor patch size 
i^/to + Tp "which then in turn determines the probability 
for binding. This line of reasoning leads to the formula 



"where a = 2td/ ^/tq, h — td is a typical diffusion 

time bet"ween binding attempts which scales ~ i/Nr, and 
Th is the homogeneous result from Eq. JSJ. Eq. ^ im- 
plies that even for vanishing receptor size the MFPT re- 
mains finite due to a finite tq. The lines in Fig. 2]sho"w 
that this equation can be fitted extremely "well to the 
simulation results. Moreover the fitted values for a and 
b agree roughly "with the predicted scaling behaviour. 

In summary, our results sho"w that the efhciency for 
initiating cell adhesion in hydrodynamic flo"w is strongly 
enhanced by increasing the number of receptor patches 
Nr, but only up to a saturation threshold. An increase of 
patch size rp leads only to a "weak enhancement of bind- 
ing efficiency. In contrast, a strong enhancement results 
from increasing the patch height rg. For example, for a 
fe"w hundred receptor patches, an elevation of ro = 10^^ 
makes initial binding already as efficient as for a homoge- 
neously covered cell. Strikingly, "white blood cells are in- 
deed characterized by such a receptor geometry, because 



they are covered "with hundreds of protrusions {microvilli, 
typical height 300 nm, corresponding to ro = 0.06) "which 
carry adhesion receptors like L-selectin at their narro"w 
tips . In general, "white blood cells operate in the limit 
of a homogeneously covered cell not only due to their 
receptor geometry, but also because during capture they 
are usually exposed to environments "with Pe « 10''— 10^ 
(a typical value for Pcz is 300) . The principle of enhanc- 
ing capture efficiency by elevation of receptor patches 
seems to be also used by other biological systems. One 
example of large medical relevance appears to be malaria- 
infected red blood cells, "which develop thousands of little 
adhesive protrusions {knobs, typical height 20 nm, corre- 
sponding to ro = 0.004) (31 . The results presented here 
do not only allo"w to understand the efficiency of cell cap- 
ture in these biological systems in a unified "way, but can 
also be used for developing corresponding applications in 
biotechnology, including adhesion-based cell or particle 
sorting. 
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